Representation Mixing and Exotic Baryons in the Skyrme Model 
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We study the effect of representation mixing in the SU{3) Skyrme model by diagonahzing exactly 
the representation-dependent part. It is observed that even without the next-to-leading order sym- 
metry breaking terms the low-lying baryon masses as well as the recently discovered Q^ and Hfo 
can be fairly well reproduced within 3% accuracy. It is also demonstrated that the mixing effect is 
not negligible in decay processes of {10}. In particular the effect of mixing with {27} is found to 
be quite large. These results are compatible with the second-order perturbation scheme. The decay 
widths are found to be sensitive to the mass values. The decay widths of {10} are estimated to be 
smaller than those of {10} by an order of magnitude due to the destructive interference between 
operators although the kinematic factors are comparable. 

PACS numbers: 12.39.Dc, 13.30.Eg, 14.20.-c 



The recent discoveries" of 8+ have generated lots of 
interesting developments in hadron spectroscopy, in par- 
ticular in understanding the exotic nature of the state. 
The state is exotic in the sense that the quantum num- 
bers cannot be explained as a system of three quarks, as 
the smallest number of quarks consistent with 0"*" is five, 
or that it cannot be classified into conventional classifi- 
cations, {8} and {10}. The lowest multiplets consistent 
with 0+ is {10} in the scheme of flavor SU{3) symmetry. 

The chiral soliton model proposed by SkyrmeQ has 
been explored theoretically andphenomenologically with 
many interesting successesQUIJIg in describing the 
properties of low-lying hadrons. The importance of the 
higher multiplets beyond octet and decuplet has been 
noticed in the chiral soliton model in treating the sym- 
metry breaking part as perturbations. The symmetry 
breaking part is not diagonal in the SU(3) multiplet 
space so that in higher order perturbation 6] 7] 8] or in 
diagonahzing the hamiltonian[9j[l0| the mass cigcnstate 
should be mixed with higher representations. For ex- 
ample, the nucleon is dominantly described by {8} but 
with non-vanishing mixing amplitudes of {10}, {27}, • • • 
and the A also has non-vanishing mixing amplitudes of 
{27},{35},.... 

The prediction 13 [iil of 9+ as the lowest state among 
the higher multiplet {10} has now been confirmed. One 
of the characteristics of 0+ as an isospin singlet and 
hypercharge 2 state with respect to representation mix- 
ing is that it has no corresponding state in the {8} and 
{27}, i.e., no representation mixing is possible. On the 
other hand more massive states in the same mulitiplet 
{10} have non- negligible mixing with other representa- 
tions and the masses and decay widths are supposed 
to depend on the mixing. The effect of mixing in sec- 
ond order perturbation has been extensively discussed 
recently [ij [i3 01 ; ^^ which the effect of mixing is found 



to be non-negligible but depends much on the parame- 
ters of the underlying effective theory. The similar ob- 
servations have been made in the exact diagonalization 
methodflal [la | for the exotic baryon masses. In this short 
note we discuss the mixing effect on the decay process 
further using the exact diagonalization method keeping 
only the chiral symmetry breaking term that is of leading- 
order in Nc- 

The effective action for the pseudoscalar mesons, which 
realizes the global SU{3)l x SU{3)ii in the Goldstone 
mode, can be written in general as 
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where ^2 and Shod are the leading kinetic term and 
the higher order derivative terms including the Skyrme 
term. Swz is the Wess-Zumino action and Ssb is an 
explicit symmetry breaking term depending on the me- 
son masses. The effective Hamiltonian after quantizing 
the "degenerate rotational mode" of the SU{2) soliton 
of hedge-hog ansatz 2] embedded into SU{3) is known to 
have the following form for Nc = 3 and B = I: 
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+ Tpr^2{SU{?,)L)~a{l~DUA)), (2) 

where C2(5C/(2)fl) and C2{SU{i)^) arc the corre- 
sponding Casimir operators [C2{SU{2)r) — J [J + 
1), C2{SU{-i)^) = \[p^ + q^ + i{p+q)+pq]. In this frame 
work we are left with four parameters, Md, /i, I2 and a, 
which should be in principle determined unambiguously 
from the effective action. In this work however we take 
them as a set of free parameters for the phenomenological 
study. 

In eq.||2I), the SU{2>) symmetric limit can be achieved 
when the last term vanishes. The mass spectrum of 
the baryon can be determined by treating the sym- 
metry breaking term in a perturbative way. One can 
also include additional terms next-to-leading order in 
chiral symmetry breaking to reproduce the low-lying 



baryon spectrum well in the first order perturbation 
calculation I 111], In this work when the haniiltonian 
is to be diagonalized, we do not include these terms 
which are of next-to leading order in the \/Nc expan- 
sion to make the analysis free from possible ambigui- 
ties due to the extra parameters in the effective the- 
ory. For the diagonalization the hamiltonian can be 
divided into two parts, representation independent (i?o) 
and dependent (_ffij) parts: 

^0 = AU + \{y y)C2{SU{2)j,) - A (3) 
Hr = ^C2{SU{Z)L)-a{l-DUA)). (4) 

Minimal extensions beyond octet and decuplet can be 
guided by considering the quark content of the baryons. 
Three quark system leads up to decuplet. With an ad- 
ditional quark-antiqauark pair for a penta-quark system, 
qqq qq, the possible representation can be extended up 
to {r0},{27},and {35}. With the constraint Yr = 1 for 
B = 1 baryon, the state vectorsJlOj for spin 1/2 baryons 
and spin3/2 baryons can be written as 
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B{J = t;)) = C^K,^) + C^oK'^) + Qrl^M .)' (5) 

■b l,T,10 \ , r^b |,T,27 \ , /-.b |,T,35 \ /q\ 
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where a(b) refer to a baryon with flavor part /^ — 
{Y,!,!^) and spin part v = (Yr, J, J3) with spin J — 
1/2(3/2). By diagonalizing the hamiltonian, Hr, in the 
form of 3 X 3 matrix for each baryon state, we can cal- 
culate the corresponding mass as an eigenvalue of the 
hamiltonian. 

The eigenvalue and the mixing amplitudes in eqs.lO) 
and (jnj are of course functions of four parameters, 
Mc/,/1,/2 and a. We fix the parameters by a best fit 
to the masses of the low-lying octet and decuplet states. 
The best fit to the mass differences can be obtained 
with the central value of/2=2.91xl0"3 MeV^and 
a = -750 MeV. Then the mass fit gives M^i = 773 MeV 
and h = 6.32 x IQ-^MeV^ It is interesting to note 
that these values are comparable to those used in the 
perturbation scheme 1Z| 13]. 

The masses in the best fit are given by 



of 0"*" and Sfo are consistent with the experimental val- 
ues within 3% accuracy. The mixing amplitudes for the 
corresponding states can be read out from the normalized 
mass eigenstates. For example, the mixing amplitudes for 
N, A and NfQ are given by 

Cf = 0.953, C^„ = 0.234, C^^ = 0.191, 
C^o = 0.877, Cfj = 0.464, C^^, = 0.125, 
Cf i"« = -0.234, C^f' = 0.970, C^r"' = 0.024, (8) 



which are comparable to those in [i3[i3Ii3- For {10}, 
it should be noted that the equal spacing rule in the first 
order perturbation is not literally respected due to the 
effects of the mixing in the 2nd order perturbation 13|. 
It is observed that there are no appreciable differences in 
the mixing amplitudes between the exact diagonalization 
scheme and 2nd order perturbation scheme, which is con- 
sistent with the high order perturbative calculations!^. 

Given the wave function in the representation space, 
eqs.© and 10, the decay width of a baryon B into a 
low-lying B' and meson (p can be obtained by evaluat- 
ing the matrix element of the baryon decay operators. 
The Yukawa coupling in general as well as the decay op- 
erator in particular, which is basically a meson-baryon- 
baryon((^i?B') coupling, has been discussed by many 
authors 18] in the context of the chiral soliton model. 
In this work, we choose an operator based on the sug- 
gestion of Adkins et al. ^j in relation to the axial current 
coupling and developed further by Blotz et al. |l3 , which 
has the form[ll|[l3: 
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where i = 1,2,3 and 6, c = 4, ...., 7. The decay amplitude 
and the decay width are given by 

Mb^b'+^ = {B'\ 6'^^^ \B), Tb^b'^ = K-A^ (10) 

where A^ = |Mp/3p^ and iiT is a kinematic factor : 



M{N) = 939, M(A) 

M(S) = 1345, M(A) 

M(Sio) = 1506, M(fi) 

M(iVi-o) = 1705, M(I]ro 



= 1108, M(S) = 1226, 

= 1231, M(Sio) == 1385, 

= 1638, Af (6+) = 1570, 

== 1811, Af(Sro) = 1818.(7) 



One can see that the masses for the low-lying octet 
and decuplet are reasonably well reproduced in the 
exact diagonalization method, with results that are 
comparable also to those obtained in the perturbation 
scheme (either in the first order p?7l| or in the second order 
perturbation|12j'). It is found that the estimated masses 
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Here Afs's(m) are the corresponding masses of the 
baryons(mesons), M s are the mean masses of the 
multiplet. We take Ms = 1154.5 MeV, Mio = 
1436 MeV, Mj^ = 1726 MeV. 

The decay amplitudes of the baryons can be calculated 
in a straightforward way and result in lengthy formulae. 
For example, the amplitudes squared for A ^ A'^-I-tt and 



0+ -^ N + K are given by 
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where Giq — Go + 9^1,627 — Go — 2Gi,G27 — Gg ^ 
2Gi,F35 = Go + GijGio = Go — Gi — ^G2,Hiq — 
Go - |Gi + iG2,ip7 = Go + iiGi + ^G2. Introduc- 
ing a parameter p [13 as Gi = p Gq, we take Gq and 
p as parameters in this phenomenological analysis. We 



find a p and Go that are consistent with the overall fit 
to the experimental values of the widths of the decu- 
plet. The overall fit is obtained with Go = 17.5 and 
p — .5. The decay width is found to be quite sensi- 
tive to the masses of the particles involved in the decay 
process. This is because the kinetic part is very sensi- 
tive to the masses. We calculate the possible range of 
the calculated widths by allowing 3% variations of the 
masses. As shown in the parenthesis in Table Q] the ki- 
netic terms K and therefore the decay widths are chang- 
ing in a relatively large range even with 3% variation with 
masses. On the other hand by allowing ±3% variation 
in masses, reasonably well reproduced in this model, one 
can explain the experimental values of {10} decay widths 
within the right range. Now given the set of parameters 
determined by the low-lying baryons, one can make the 
prediction for the decay widths of exotic {10} baryons. 
In this work we adopt the parametrization for G2 as in 

E 

widths are given in Table HI] 



G2 — ( 3F/D+5 ) ^P ~^ 2)Go. The estimated decay 



J 



TABLE I: {10} ^ {8} + <p 



Decay 



K" 



A^ 



T^a 



A^N + TT 


0.33 (0.13~0.64) 


Sio ^ A + TT 


0.17 (0.04-0.44) 


El„ ^ E + TT 


0.001 ( < 0.18) 


Hio ^ H -1- TT 


0.01 ( < 0.22) 



367 
177 
43 
135 



121 (46 ~ 233) 
31 (8 ~ 79) 
0.70 ( < 7.9) 
1.2 ( < 30) 



115 - 125 
34.7 
4.73 
9.9 



"Values in the parenthesis are obtained with ±3% mass variations. 



Decay 



TABLE II: {10} -^ {8} + ip 



K" 



42 

^(bestfit) 



^(Gi=0,G2=0) 



42 
^{without{27}) 



e+ ^ N + K 


0.52 (0.15 ~ 1.04) 


7.60 


165 


29.00 


4 (1.2 ~ 7.9) 


Efo -^ H -1- TT 


0.66 (0.42 ~ 1.3) 


39 


125 


16 


26 (16 ~ 50) 


Sio ^ E + if 


0.23 (0.06 ~ 0.86) 


17 


41 


17 


4 (0.97 ~ 14) 


Nio ^N + 7v 


3.3 (2.5 - 4.2) 


0.61 


15.2 


14 


2 (1.5 ~ 2.6) 


Nfo ->N + rj 


1.1 (0.56 ~ 1.8) 


6.6 


24 


14 


7.3 (3.6 ~ 12) 


Nfo ^A + K 


0.28 (0.01 ~ 0.67) 


3.6 


31 


0.7 


1.0 (0.03 ~ 2.40) 


Nro -^T. + K 


- ( < 0.27) 


0.40 


45 


2.3 


- ( < 0.11) 


Efo -^N + K 


2.4 (1.6 ~ 3.3) 


0.50 


22 


0.02 


1.2 (0.81 ~ 1.7) 


Ero ^ E -H TT 


1.3 (1.0 ~ 2.2) 


1.3 


13 


12 


1.7 (1.3 ~ 2.9) 


Efo -> E 4- ry 


0.06 ( < 0.57) 


18 


46 


34 


1.0 ( < 10) 


Ero ^ A -h 77 


0.57 (0.13 ~ 1.1) 


1.1 


22 


1.1 


0.61 (0.14 ~ 1.2) 


Efo ^S. + K 


- ( < 0.20) 


33 


70 


91 


- ( < 6.60) 


"Values in the parenthesis are obtained with ±3% mass 


variations. 









Compared to the decuplet, the decay amplitudes for 
the antidccuplct arc found to be much smaller by an or- 
der of magnitude whereas the kinetic terms are compa- 



rable to each other. It has been understood that this is 
mainly due to the destructive interference between the 
operators im. In the fourth column, the amplitudes with 



Gi = G2 = are shown, which clearly shows that the 
effect of interferences are substantially large. To see the 
effect of representation mixing particularly with {27} , 
the results without {27} mixing are shown in the fifth 
column. The overall tendency is that the non-vanishing 
mixing with {27} reduces the amplitudes [l^. However 
for the processes Sfj 



s + TT and 



^10. 



N + K, 



the mixing enhances the decay amplitudes [13, whereas 
Sfo — > S + ii" and SfQ -^ A + 77 are found to be insensi- 
tive to {27} mixing. The values in parenthesis are those 
with ±3 % variations of the baryon masses. According 
to our calculated masses, the process A^fg ^' T, + K and 
E]7o -^ E + K are beyond the threshold in the best fit. 

In this work, we discussed the effect of representation 
mixing obtained in SU{S) Skyrme model by diagonalizing 
the representation dependent part in the hamiltonian re- 
sulting from quantizing the rotational mode. It is shown 
that even without the ncxt-to- leading order (in Nq) sym- 
metry breaking terms the low-lying baryon masses can 
be fairly well reproduced by allowing the mixing with 
higher representation. One of the major differences in 
the mass results obtained in the exact diagonalization 
method compared to the first order estimation jljj is that 
there is a deviation from the equal spacing rule with hy- 
percharge in the {10} multiplet . It is due to the non- 



negligible mixing with other representations^!- It is also 
observed that the mixing effect is not negligible in the 
decay widths. The effect of mixing with {27} is found 
to be particularly large. These results are consistent 
with the second-order perturbation scheme, where hi ghe r 
order corrections are found to be relatively large ^^■ 
Although the decay-width estimations in this work are 
based on a specific form of the decay operator 4] 19], the 
observation that the results of the exact diagonalization 
method and the second-order perturbation scheme are 
consistent with each other demonstrates that the higher 
order corrections beyond the second order might not be 
important in numerical estimations. However, it should 
be noted that the exact diagonalization can make more 
sense only when the hamiltonian to be diagonalized is as 
complete as possible at least for the symmetry breaking 
part. 
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